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The Higgs boson H has the largest coupling to the top quark t among the standard model (SM) 
fermions. This is one of the ideal places to investigate new physics beyond SM. In this work, we study 
the potential of determining Higgs boson CP properties at the LHC and future 33 TeV and 100 TeV 
pp colliders by analysing various operators formed from final states variables in tiH production. The 
discrimination power from SM coupling is obtained with Higgs boson reconstructed from H —>■ 77 
and PI —>■ bb. We find that tfbb process can provide more than Sa discrimination power with 300 
fb~^ integrated luminosity in a wide range of allowed Higgs to top couplings for the LHC, the 33 
TeV and 100 TeV colliders. For tf~f~f the discrimination power will be below 3cr at the LHC, while 
for 33 TeV and 100 TeV colliders, more than 3cr sensitivity can be reached. 


I. INTRODUCTION 

Since the discovery of the Higgs boson at the LHG by 
the ATLAS and GMS Gollaborations [l|, , experimental 

measurements have been focusing on the determination 
of its consistency with standard model (SM) expectation. 
Within the current experimental and theoretical uncer¬ 
tainties, the discovered Higgs boson is in agreement with 
the predictions for a SM Higgs boson considering its spin 
zero property ii- The measurement of Higgs coupling 
to fermions through 66 [^,Q, T+r“ @,[1], ttiT [g.[^ and to 
gauge bosons through WW [IS[II1 j [i3,[i 3 have also 
been performed with current LHG data and will be mea¬ 
sured more precisely in the future run. The measurement 
of Higgs self-coupling is also necessary to reconstruct the 
scalar potential of the Higgs doublet field, while due to 
the smallness of the signal and large QCD backgrounds 
the probing of the Higgs self-coupling has to await for a 
high-luminosity LHG (HL-LHC), 33 TeV, or 100 TeV pp 
collider Urn- 

It is important to stress that so far all experimental 
determinations of the Higgs CP properties have been ob¬ 
tained from the Higgs to vector boson couplings with lep¬ 
ton final states. The hypothesis that Higgs boson being a 
pure pseudoscalar state has been excluded by the present 
data since the 2i2iy and 4£ signals from H —>■ WW 
[sini and ZZ [H, [T^ decays have been observed. In 
particular, the angular distributions of the lepton pairs 
in the H ^ ZZ channel are sensitive to the spin-parity 
of the Higgs boson. Although the present LHG data 
strongly prefer that iL is a = 0+ state, it is not yet 
excluded that the Higgs boson has a pseudoscalar compo¬ 
nent SHI- Using the present measured signal strengths 
in various Higgs search channels, theoretical studies on 


the Higgs CP properties have been performed, and con¬ 
straints on the mixing angle of the GP-even and GP-odd 
component are given model dependently [itI - I^ . There 
are also investigations on how and to what extent to pin 
down the Higgs CP mixing angle via other Hi ggs decay 
channels at the LHG and future linear collider [27H^ . 

Among all the Higgs production channels, the investi¬ 
gation of the tiH process plays a complementary role in 
Higgs characterization in the sense that it is very sensi¬ 
tive to the relative magnitudes of the GP-even and CP- 
odd top-Higgs Yukawa coupling coefhcients. Although 
the tiH production cross section is only around 1/200 
of the inclusive Higgs production for with a mass of 
125 GeV at run I LHG and suffers from large multiplic¬ 
ity of the objects in the final state from top pair decays, 
the HL-LHG and proposed future hadron colliders with 
higher collision energy, the high energy LHC (HE-LHC), 
SppC, etc., are capable of enhancing the tiH event rates. 
In this paper, we will revisit certain operators defined 
from the momenta of the tiH process, which are sensi¬ 
tive to the coefficients of different GP-mixed states. We 
choose two particular Higgs decay modes 77 and 66 for 
clearness of signal and the wealth of event rates respec¬ 
tively, and study the behaviour of weighted moments at 
the LHC and the future hadron colliders. We find that 
these operators for the corresponding 6677 and tibb pro¬ 
cesses have a good discrimination power for a CP-even 
Higgs from a GP-mixed state. 


II. EFFECTIVE COUPLINGS FOR THE 
GP-MIXED HIGGS STATE 

To accommodate possible deviations of Higgs to top 
quark tiH coupling in both strength and CP properties, 
we parametrize the tiH coupling in the following form 

Cnti — —t(cos^ -\- C^sinffjtH, 
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where v is the vacuum expectation value (vev) of Higgs 
field, and sin ^ = 0 and cos ^ = 0 corresponds to the CP- 
even and CP-odd coupling coefficient respectively. In 
particular, the SM Higgs boson has a = 1 and cos^ = 1. 
The angle ^ is an indication of the degree of mixture CP 
in the physical Higgs H. This form of Higgs fermion 
interaction will be assumed also for other fermions in our 
discussions. 

The dominant Higgs couplings to vector gauge bosons, 
W and Z will come from the CP-even component of H 
in the following form 

Chvv = a cose . (2) 

The main production channel of Higgs boson at the 
LHC is via gluon fusion process. The effective vertices 
for a CP-mixed Higgs state interacting with gluons at 
one loop order are given by (3ll - l35j| 



^Hgg 




H. 


(3) 


with/9 = acoseEi=6.t^a(Ti), 4® = asineEi=M-^f'('^*) 
defined as the scalar and pseudoscalar form factors for 
Hgg which retain the dominant contributions from top 
and bottom quarks, and = mfj/Amf. G^^ is the 
gluon field strength and G^j, is the dual of given by 
(f/2)e^j/a/3G“^. The two form factors Fair) and Pb(T) 
can be expressed in terms of scaling function /(r) as 


Fa{T)=T l[l-h(l-T ^)/(t)] , Pb(r)=T ^ f {t) ,{F) 
where 

if T< 1 , 


/(t) = 


1 [sin 


1-1 



n 2 


(5) 


The Higgs to diphoton effective coupling is formulated 


as 




ji p p^^ -I- P F 




( 6 ) 


where F^^ and P^j/ are the photon field strength and 
its dual, respectively. The form factors /Jj, are given 

by 12 = a cos ^[2 J2i^b,t^cQFa{Ti) - Fi{t^^)], I2 = 
2asin^ J2i=b t ^cQi^bin) with Nq = 3 the colour factor 
for quarks and Qi the electric charge of the corresponding 
particle. Fi{t-^) is from W boson loop contribution with 

Pi(r) = 2 + 3r-i + 3r-42 - (7) 


III. PRODUCTION OF ttH AT pp COLLIDERS 

Deviations of Higgs boson coupling to top quark from 
SM one can show up in different ways. We study how 
to identify such deviations using pp —>■ tiFIX with the 


FIG. 1: The cross section ratio ct/ctsm versus CP mixing 
angles for t/yy and ttbb processes at ^/s= 13 TeV, 14 TeV, 33 
TeV, and 100 TeV. In the figure, the parameter a is set to be 

1 . 


Higgs boson identified by // —>■ 66 and F[ ^ since the 
former has the largest branching ratio and the latter has 
clearer signals compared with other decay modes. 

Both ATLAS [§, and CMS @ have performed ex¬ 
perimental searches on the tiFl production process, with 
H decays into yy and 66 separately. For the ttbb process, 
collecting the i/s = 8 TeV data final state are categorized 
according to their jet and 6-tagged jet multiplicities. The 
observed cross section upper bound is d.lcrsM [1^, and 
3.3 (tsm 0- For the ttyy process, the observed upper 
limit on cross section is 6.7 (tsm 0- We still need to wait 
for more accurate data to draw conclusions about the 
production cross section. 

The interaction in Eq.(IT|) can cause deviation in ttH 
cross section from SM prediction and in the Higgs decays. 
We have evaluated the ttH cross sections using Mad- 
Graph5_aMC particularly heft model where the Higgs 
effective theory is included [sj with a mild rapidity cuts 
on the top and antitop quarks |77 j yj < 4. The results 
with the parameter a = 1 and range of ^ from 0 to tt 
and collision energy ^/s = 13 TeV, 14 TeV, 33 TeV and 
100 TeV with = 125 GeV are shown in Fig. [T] Here 
we would like to point out that one should consistently 
take the high-order correction to the cross section into ac¬ 
count. For ttH process, the NLO prediction for a Higgs 
boson mass of 125 GeV has cross section around 86 fb, 
130 fb, 611 fb and 33700 fb [H-Iiil at = 7 TeV, 8 
TeV, 14 TeV and 100 TeV. Estimated using CTEQ6.6, 
the NLO corrections are positive. This is higher than 
that obtained with the cross section of about 72 fb, 107 
fb, 476 fb and 24010 fb using the LO results. A //-factor 
of 1.19 - 1.40 [11,0,143 should be applied when we use 
the values shown in Fig. [TJ 
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With a fixed at 1, one obtains the cross sections for 
the SM for ^ = 0. When going away from ^ = 0, one 
can clearly see that the cross sections become consider- 
ablely smaller than those of the SM predictions. One 
wonders this fact can be used to distinguish the SM UH 
coupling with that from a beyond SM. Global fit on the 
CP-violating phase using the observed rates oi gg ^ H 
and —>■ 77 with the LHC 7 and 8 TeV data can be 
translated into constraints on ^ as |^| < O.TStt at 99.73% 
C.L assuming a = 1 [l^. It has been shown that if the 
LHC measurement of aitiH) is with an accuracy of 20%, 
it is indicated that |f| < 0.177r can be obtained [^ . 
When the cross section is measured to good precision, 
one can obtain some important information about the 
interaction. 

If experiments obtain a cross section different from that 
with a = I and ^ = 0, one immdediately knows that 
beyond SM physics is needed. However, should a cross 
section close to the SM will be obtained, one cannot rule 
out the possibility of beyond the SM ttH coupling given 
in Eq. ©• This is because that even with a ^ nonzero and 
therefore smaller cross sections as shown in Fig. [T] if the 
parameter can be varied from I, one can still adjust the 
value of to be inverse of ct/ctsm to obtain similar val¬ 
ues as SM ones. One, however, note that just measuring 
cross sections, one will not be able to have information 
about Higgs CP properties. Therefore it is desirable to 
find ways to distinguish SM from beyond SM ttH cou¬ 
plings independent of the overall scaling parameter a and 
also to provide information about Higgs CP properties. 
In the rest of this paper, we show that certain weighted 
moments in pp ttHX process are sensitive to the rela¬ 
tive magnitudes of the CP-even and CP-odd interactions 
in tiH. Defined from the ratios of integrated operators 
to total cross sections, these weighted moments are also 
partially free from NLO corrections, where the K factors 
can be factored out. 


IV. WEIGHTED MOMENTS IN UH PROCESS 


To achieve this, we take the same approach as those 
proposed in Ref. to study operators in the following 
formed from final product variables which are sensitive 
to cos^ ^ — sin^ f 


(^1 

(^3 

(^5 


{pt X n) ■ {pt X fi) 

|(pj xn)-(KXn)r " 

{pt X n) ■ (pt X n) 

rp rp , L /4 - 

Pt Pf 

PtPi .n - pM 

pfpr "“iptiipti’ 


. PtPi 
■ \pfp^\ 

{pt X fi) ■ {pi X n) 

IPtllPfl 


( 8 ) 


where pjij. denote the magnitudes of the t and t transverse 
momenta, n is a unit vector in the direction of the beam 
line and defines the 2 : axis, while x axis is chosen to be 
any fixed direction perpendicular to the beam. There 


may be other operators but we will restrict ourselves to 
the six operators in Eq. as examples. 

As mentioned before that the unknown parameter a 
may make interpretation of information extracted diffi¬ 
cult, it is therefore desirable to find observables that are 
independent of the parameter a. To this end we define 
the following observables by taking ratios to remove the 
a dependence, 

^ J[0,]{da{pp -» tiXX)/dR}dR 

J{da{pp^tiXX)/dR}dR ’ 

where XX represents Higgs decay products and R is the 
phase space of tiXX process. The values of ag are shown 
in the upper panel of Fig. [2] One can see that the changes 
of as against ^ is significant providing hope to distinguish 
Higgs boson coupling with top quark with mixed CP 
component. 


V. DISCRIMINATION POWER FOR ttyy AND 
ttbb PROCESS 


We now study to what accuracy for a given it can 
be measured. One must consider errors associated with 
the variable when both signal and SM background are 
taken into account. The statistic error for is given by 


(5q!c — 




/3s - - “iasas P 


1/2 


,( 10 ) 


with S and B the total number of events for signal and 
background process respectively. The variable PsH for a 
given operator Ot is defined as 


1 _ lPi?{d^iPP ttXX)/dR}dR 
* J{da{pp ^ ttXX)/dR}dR ’ 

and shown in the lower panel of Fig l2] for different oper¬ 
ators varing the CP mixing angle ^ at different collision 
energies. Note that l3g{Oi,2) are equal to 1. 

To quantify the ability to distinguish the SM pure 
scalar case from any CP-mixed Higgs state of each oper¬ 
ator, the discrimination power D is defined 


|a|^(^ = 0)-4(0| 

= 0 ) 


( 12 ) 


We vary the mixing angle ^ from 0 to tt when a series of 
discrimination powers can be obtained accordingly from 

Eq.dUD. 

For a full consideration, we must include both signal 
and background for Higgs and top pair decay. Since H —>■ 
bb has the largest branching ratio and iL —^ 77 has clearer 
signal, we will study the behaviours of ag for the ttbb and 
tt 77 processes aiming for the potential of these operators 
to discriminate a nonzero ^ compared with SM coupling 
in light of new data available at different energies for a 
pp collider. 
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FIG. 2: The as and Ps as functions of ^ for pp — >■ ttHX 
process at ■\/s= 13 TeV, 14 TeV, 33 TeV and 100 TeV. 


We compute the cross section of the ^77 and ttbb pro¬ 
cesses u sing MadGraph5_aMC for both signal and back¬ 
ground Cross sections for ttjj production are ob¬ 
tained after basic cuts \rit,t,'y\ < 4, > 25 GeV and the 

mass window cut \Mry^ — < 5 GeV, with = 125 

GeV. For the ttbb process, we impose the basic cuts with 
\%,iAb\ < 4 and > 25 GeV. 

If we further include the ti pair decay, a total factor 
of the change on the cross section cr and discrimination 
power D can be estimated by the corresponding decay 
branching ratios. We choose top quark decays semilep- 
tonically and antitop quark decays hadronically in order 
for the identification of t and t and also the reconstruction 
of the transverse momentum of t and t. Therefore, the 
decay products from top pair will be 26-jet + 2 j 


the reduced signal cross section can be estimated by n = 
B{t ^ bWfB{W ^ lv)B{W ^ jj)cs x i x |, with 
n ~ 0.18 for 6677 process. While for ttbb cross section, 
an additional suppression factor from 6 -tagging efficiency 
60%) should be included leading to n{tibb) ~ 0.0648. 
Gonsequently, by including top pair decay the value of 
D will be affected by a total suppression factor of -Jn ~ 
0.42 (0.25) for 6677 [ttbb) from Ea. lfT^ . 

Given a^, /3g, a^, fig and background event number 
B, it is straightforward to compute the required signal 
event number S for a SM-like Higgs to achieve certain 
value of D. For example, if one wants to discriminate 
ttH coupling of ^ = 0 from p = O.Stt at the 3a statisti¬ 
cal level, the prediction from 6677 and ttbb are different. 
The 6677 background event number is H ^ 50 (2280) at 
= 13 (100) TeV, estimated from the operator with 
the largest discrimination power O 3 we find that D = 3 
requires signal event number to be around 750 (1950). 
The simulated value for LO signal rates are about 160 
(9870). Therefore, using ttjj channel it is difficult to 
discriminate the ttH of ^ = 0 coupling from a ^ = O.Stt 
mixed coupling at 3a statistical level at 13 TeV. How¬ 
ever, a future pp collider at 100 TeV can produce enough 
signal events to achieve the discrimination at larger than 
3a level. The situation is better in ttbb channel. The SM 
ttbb background has a event number of around 8.0 x 10 ® 
(7.8x10^) at ^/s = 13 (100) TeV, estimated from the 
operator with largest discrimination power Oq we find 
that D = 3 requests signal events number to be 4.1 x 10"^ 
(4.3 X 10®). The real LO signal event number is around 
7.1 X 10^ (4.4 X 10®), which is significantly lager than re¬ 
quired. Gonsequently, from the ttbb process, one is able to 
discriminate the tiH coupling with ^ = 0 from a ^ = O.Stt 
CP- mixed state at more than 3a statistical level with any 
collision energy larger than 13 TeV at pp collider. 

The figures [3] and [4] show the corresponding D values 
for ^77 and ttbb final state after including the above de¬ 
cay and 6 -tagging efficiencies with the CP mixing angle 
running from 0 to tt with collision energy at 13 TeV, 
14 TeV, 33 TeV and 100 TeV. For both processes, the 
number of events are obtained with integrated luminos¬ 
ity assumed to be £ = 300 fb~^ and an enhancing factor 
of -x/lO ^ 3.2 on D will appear if we consider the HL- 
LHC with C =3 ab“^. Gomparison between each two 
operators shows that Ci,C 3 ,C 4 , and Oq are most use¬ 
ful balancing the sensitivity and errors, among which O 3 
probably performs the best. O 3 is also easy to be con¬ 
structed which is irrelevant to the z-component of top 
and antitop momentum 

The results shown in Fig. [3] tell that for 6677 process 
at \/s= 13 TeV and 14 TeV the discrimination powers 
for operators Oi-Oq are generally less than three in the 
allowed CP mixing region. When the integrated lumi¬ 
nosity is increased to be 3 ab“^ at HL-LHC, D > 3 can 
be reached in the region |^| > 0.3Tr for all operators. It 
is much more optimistic at ^/s = 33 TeV and 100 TeV, 
which shows that all the operators can give large enough 
statistical significance in the region |^| > 0.2Tr. The oper- 
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FIG. 3: The discrimination power D versus CP mixing angles ^ for each operator for tt'y'y process at v^=13 TeV, 14 TeV, 33 
TeV and 100 TeV . 


ator O 3 gives the best discrimination power irrelevant of 
the collision energy and CP mixing angle. All the oper¬ 
ators reach maximal discrimination power at |^| = O.Stt. 


Comparing the ttbb process with ttyy, we want to em¬ 
phasize, as can be seen from Fig. 01 that H ^ bb de¬ 
cay modes shows better discrimination powers due to the 
large event rate of the Higgs to bb decay branching ra¬ 
tio. For tibb process, the operator O 3 and Oq give the 
best discrimination powers. Note that if we make full 
analysis considering final state jets and leptons from top 
quark decay, and impose further cuts on the jet multi¬ 
plicities, the value of discrimination power D will appear 
even better than present numbers. 


VI. SUMMARY AND DISCUSSIONS 

The Higgs boson H has the largest coupling to the 
top quark t. The ttH interaction can be sensitive to the 
investigation of new physics beyond SM. In this work, we 
have studied the potential of determining Higgs boson 
CP properties at the LHC and future 33 TeV and 100 
TeV pp coliders by analysing various operators formed 
from final states variables in ttH production. 

A CP violating Higgs coupling to top quark can cause 
significant deviations from SM predictions for pp —>■ 
tiHX. We have evaluated the tiH cross sections us¬ 
ing MadGraph5_aMC particularly heft model. We find 
that the cross sections become considerably smaller than 
those of the SM predictions when CP violating compo- 
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FIG. 4: The discrimination power D versus CP mixing angles ^ for each operator for ttbb process at \/s = 13 TeV, 14 TeV, 33 
TeV and 100 TeV. 


nent gets bigger. This fact may be used to distinguish 
the SM ttH coupling with that from a beyond SM. But 
this may be caused by an overall scaling factor to the 
coupling. To further identify the relative magnitudes of 
the CP-even and CP-odd interactions in ttH, certain 
sensitive weighted moments in pp —^ ttHX process can 
help to achieve this. We obtained results for discrimi¬ 
nation power for several operators with the Higgs boson 
identified by P —?► 66 and H ^ For an integrated lu¬ 
minosity is 300 fb~^, we find that, with ttH —>■ 6677 , the 
discrimination power will be below Scr at the LHC, while 
for 33 TeV and 100 TeV colliders, more than 3a sensi¬ 


tivity can be reached. On the other hand, ttH —>■ tibb 
process can provide more than 3cr discrimination power 
in a wide range of allowed Higgs to top couplings for the 
LHC, the 33 TeV and 100 TeV colliders. 
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